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OPTICAL  GUIDING  IN  THE  SEPARABLE  BEAM  LIMIT 


T.  M.  Antonsen,  Jr.  and  B.  Levush 
Laboratory  for  Plasma  and  Fusion  Energy  Studies 
University  of  Maryland 

ABSTRACT 

The  nonlinear  theory  of  optical  guiding  in  an  FEL  Amplifier  is 
developed  for  the  case  in  which  the  spatial  dependence  of  the  current 
source  term  in  the  wave  equation  can  be  separated  into  the  product  of  a 
function  of  radius  and  a  function  of  axial  distance.  Such  a  separation 
can  be  motivated  if  either  the  betatron  wavelength  is  shorter  than  other 
lengths  of  interest  (synchrotron  wavelength,  vacuum  Rayleigh  length)  or 
if  the  radiation  waist  exceeds  the  beam  radius.  In  this  limit  with  the 
choice  of  a  Gaussian  profile  for  the  electron  beam  density  the  wave 
equation  can  be  solved  exactly  and  the  radiation  field  felt  by  the 
particles  can  be  expressed  as  a  one-dimensional  convolution  of  the 
current  source. 

With  the  given  expression  for  the  radiation  field,  the  equations  of 
motion  can  be  solved  in  the  trapped  particle  regime.  Requiring 
consistency  between  the  particle  motion  and  the  fields  yields  expressions 
describing  nonlinear  guided  states.  The  adiabatic  evolution  of  these 
guided  states  in  the  presence  of  a  tapered  wiggler  is  determined  by 
conservation  of  the  electrons'  action  and  total  (field  +  electron  beam) 
energy.  Using  these  relations  the  growth  of  the  radiation  waist  as  the 
beam  is  decelerated  can  be  calculated. 

I.  INTRODUCTION 

Optical  guiding  in  free  electron  lasers  can  become  a  complicated 
numerical  problem  in  the  nonlinear  regime  due  to  the  three  dimensional 
structure  of  the  wave  fields.  Analytic  treatments  for  the  linear 

regime  of  exponential  gain  have  been  developed.  For  nonlinear 

simulations  the  representation  of  the  field  can  be  accomplished  either  by 
finite  difference  methods  or  expansions  in  a  set  of  basis  functions. 
Invariably  questions  are  raised  concerning  the  accuracy  of  these 
representations.  The  purpose  of  this  paper  is  to  show  that  with  a  few 
simple  heuristic  approximations  to  the  current  source  term  in  the  wave 


equation  the  wave  equation  can  be  solved  exactly.  The  result  is  that  the 
particle  motion  can  be  treated  as  in  a  one  dimensional  model  and  many 
analytic  results  describing  the  effects  of  diffraction  can  be  obtained. 
Furthermore,  this  model,  whose  solution  can  be  obtained  analytically,  can 
serve  as  a  test  bed  for  the  other  methods.  The  model  we  present  is 
similar  to  one  suggested  by  Moore  in  the  strong  diffraction  limit 

We  begin  with  a  consideration  of  the  wave  equation  for  the  signal 
vector  potential  Ag  which  is  appropriate  for  amplifiers, 
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where  6j(x^,z)  gives  the  complex  amplitude  of  the  component  of  the  beam 
current  density  oscillating  in  phase  with  the  signal  wave.  The  current 
density  will  be  localized  in  the  transverse  dimension  to  an  area  the  size 
of  the  electron  beam.  Furthermore,  we  can  expect  if  all  the  electrons 
are  interacting  strongly  with  the  signal  field  the  z  dependence  of  6j  at 

Accordingly , 

we  will  make  the  approximation  that  the  axial  and  transverse  dependence 
of  6j  can  be  separated. 


different  radial  positions  x^  will  be  strongly  correlated 
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where 
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me  y 


is  the  jitter  velocity  induced  by  the  wiggler  field,  vz  and  I  are  the 
speed  and  current  of  the  electron  beam,  and  we  have  taken  the  electron 
beam  to  have  a  Gaussian  profile  in  radius  characterized  by  a  width  r^. 
Here  we  have  assumed  a  planar  wiggler.  In  Eq.  (2)  the  function  s(z) 
gives  the  axial  dependence  of  the  current  density  and  will  be  determined 
by  the  dynamics  of  the  interaction  of  the  beam  particles  with  the 
ponderomotive  wave.  By  analogy  with  the  one  dimensional  problem  we  write 


s(z) 
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where  Mz)  represents  the  phase  of  a  particle  in  the  beat  wave  potential, 
and  the  angular  average  Is  over  initial  particle  phases. 


This  separation  is  probably  not  justifiable  rigorously  except  when 
the  radiation  waist  is  much  larger  the  electron  beam  radius  and  details 
of  the  radial  dependence  of  the  current  source  are  washed  out.  This  is 
essentially  the  argument  advanced  by  Moore. ^ 

The  particles  in  our  model  evolve  according  to  the  standard  one¬ 
dimensional  pendulum  equation. 


8 


.  .  o>v  ,  qA 

—  (5v  +  v  )  =  -  L_  r—2)  I  r_L 
dz  v  T  TR'  lv  cJ  2  L  2 

z  me 


3ty 

d  ,  r  3  , 

d?  ♦  '  -  f  s?7  ( 


s  i^ 
e 


_* 

qA, 


s 


me 


Vz(Vz) 


(4a) 

(4b) 


where  yR(z)  is  the  so  called  resonant  value  of  y  that  is  required  to  keep 
particles  in  phase  with  beat  wave  in  the  presence  of  tapered  wiggler 
parameters,  and  6y  is  the  deviation  of  a  particle's  actual  y  from  the 
resonant  value.  The  quantity  Ag  is  the  effective  field  that  a  beam 
particle  feels  and  represents  some  radial  average  of  the  actual  field 
Ag(r,z).  In  our  separable  beam  approximation  there  is  no  "first 
principles"  prescription  for  making  this  average.  However,  it  can  be 
verified  afterward  that  if  one  makes  the  choice 
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then  a  conservation  of  energy  relation  exists  for  the  system. 

Accordingly,  we  pick  Eq.  (5)  to  define  A  (z). 

s 

The  system  of  equations  (4)  can  be  put  into  the  dimensionless  form 
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where  the  normalizations  of  5i|>,  and  a  are  given  in  Ref. 
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by  the  Rayleigh  length  LR  =  r^w/c  anc* 
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9  with  L  replaced 


(7) 


measures  the  deceleration  due  to  tapering.  The  wave  field  satisfies  the 
Schrodinger  like  equation, 
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measures  the  current  and  1^  =  mcJ/q  =  1.7  x  10  Amps  is  related  to  the 
Alfv€n  current,  and  P  =  r/rg  is  the  dimensionless  radius.  The  only  other 
parameters  entering  the  system  are  the  dimensionless  beam  current  I,  and 
the  distribution  of  injected  momenta  f(p,5=0). 

From  this  system  of  equations  we  may  derive  the  following  energy 


conservation  law 
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II.  SOLUTION  OF  THE  WAVE  EQUATION 


We  now  turn  to  the  wave  equation  (8).  Because  the  radial  dependence 
of  the  source  is  prescribed  to  be  Gaussian  we  can  solve  for  the  wave 


field  in  terras  of  a  Green's  function  G(p ,£  -£')  which  satisfies 


V^G  +  2i  P-).  6(g  -  5')  . 


The  resulting  expression  for  the  wave  field  is  then  expressed  as  a 


convolution, 


a(p,5)  -  I  /  d£ '  G(p,5  -  S')S(5')  +  a^(p,5) 
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where  a^(p,£)  represents  the  homogeneous  vacuum  solution  corresponding  to 


the  initial  (injected)  wave  field.  Equation  (11)  is  easily  solved  to 


yield. 
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and  the  radial  average  in  Eq.  (5)  can  be  performed  to  obtain  a(£) 
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Thus,  we  see  that  the  effective  field  felt  by  the  particles  can  be 


written  as  a  one-dimensional  convolution  of  the  current  source. 
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The  kernal,  (l  +  i(£  -  £'))  ,  has  an  obvious  physical  interpreta¬ 
tion.  For  small  values  of  (£  -  5')  the  wave  field  is  just  an  integral  of 
source.  This  is  the  result  that  would  be  obtained  if  we  had  dropped  the 
transverse  derivatives  in  the  wave  equation  and  represents  the  fact  that 
nearby  contributions  of  the  source  to  the  field  are  the  same  as  if  there 


were  no  diffraction.  For  large  values  of  5  -  (recall  distance  is 

normalized  to  the  Rayleigh  length)  the  Kernal  decays  weakly  indicating 

the  fact  that  distant  contributions  of  the  source  to  field  are  reduced  by 

diffraction.  We  note  that  this  decay  is  very  weak  and  that  replacing  the 

transverse  derivatives  in  Eq.  (7)  by  a  constant  loss  rate  which  yield  an 

exponential  decay  of  the  kernal  greatly  overestimates  the  decay. 

It  is  interesting  to  examine  the  amplitude  of  the  field  that  is 

predicted  when  the  source  has  an  exponential  dependence  on  axial 

iv£ 

position,  S(£)  =  S^e  .  Such  a  dependence  certainly  occurs  in  linear 
theory.  We  shall  see  also  that  this  dependence  occurs  for  steady 
nonlinear  states  as  well.  It  is  found  that  as  (  +  » 
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Using  Eq.  (6b)  we  can  calculate  the  rate  at  which  energy  is 
extracted  from  the  particles  to  maintain  this  field, 

I<p>  =■  -  [aS*  -  a*S] 

H  2  (17) 

-  -  i  T  |S0|2  el(v  -V**[E(v)  -  ( E(v ))  *]  . 

Thus,  the  imaginary  part  of  E(v)  gives  the  rate  at  which  energy  is 
extracted  from  the  particles.  For  v  real  and  positive  E(v)  is  also  real 
and  the  field  is  maintained  in  steady  state  without  energy  being  supplied 
by  the  particles.  If  v  is  real  and  negative  we  can  manipulate  the 
contour  in  Eq.  (16)  to  show 

4n  j  Q  1  -  y  4 

where  P  indicates  the  principal  part  of  the  integral.  The  particles  lose 
energy  according  to 
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This  energy  is  radiated  away  from  the  beam  and  represents  a  situation  in 

which  optical  guiding  is  not  present. 

The  preceding  discussion  shows  that  the  key  requirement  for  optical 

1  8 

guiding  in  our  model  is  that  v  >  0.  ’  Here  our  definition  of  an 

optically  guided  state  is  one  that  can  be  maintained  without  the 

continuous  transfer  of  energy  from  the  particles  to  the  fields. 

Physically,  the  condition  v  >  0  is  simply  the  requirement  that  the  axial 

phase  speed  of  the  perturbed  current  source  be  less  than  the  speed  of 

light.  In  this  case  the  fields  outside  the  electron  beam  are  evanescent 

and  field  energy  is  confined  to  the  vicinity  of  the  electron  beam. 

Equation  (17)  may  also  be  used  to  calculate  the  field  energy  stored 

in  an  optically  guided  state  with  v  >  0.  We  imagine  such  a  state  is 

created  adiabatically  and  let  v  =  v  +  iv,  with  v,  <  0  and  |v  I  »  |v  I. 

r  i  i  1  r 1  1  i ' 

We  find  by  using  Eqs.  (9)  and  (17) 
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Integrating  from  £ 


*  -  00  to  5  and  letting  v  approach  zero  we  have 
“  o  ,  3  E/3v 

/  2tt p  dp  |  a  |  =  4 - ^ —  |  aQ  |  2  • 


(18) 

0  "  E* 

Thus,  we  can  now  express  the  field  energy  in  a  guided  state  in  terms  of 

the  amplitude  of  the  effective  field  aQ. 

Finally,  we  define  an  expression  for  the  beam  to  radiation  filling 

factor,  F. 

2,  ,Z 
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or  for  steady  states, 
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The 


Using  Shwartz's  inequality  one  can  show  from  Eq.  (19a)  F  <  (2u) 
limit  F  =*  (2ir)  is  achieved  in  steady  state  for  large  v  which  follows 
from  the  asymptotic  expansion  of  r-q.  (16).  This  conforms  with  our  inter¬ 
pretation  of  the  kernal.  For  large  v  contributions  to  the  source  from 
nearby  positions  dominate.  As  a  result,  diffraction  is  not  important, 
and  the  beam  and  radiation  are  perfectly  matched.  On  the  other  hand,  for 

-1  2 

small  v  we  obtain  F  ~  (2x)  v£n  v.  In  this  case,  due  to  diffraction  the 
electron  beam  occupies  a  much  smaller  area  than  the  radiation  field. 
Note,  however,  that  F  decreases  slowly  with  v.  For  v  =  .01  we  find 


J  l.l  F.n.ri.ri.nl  i.n.l  i.l 
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2tt F  =  0.27.  This  points  out  that  diffraction  is  a  weak  effect. 

III.  LINEAR  THEORY 


In  this  section  we  will  analyze  the  linear  regime  described  by  our 
separable  beam  model.  To  aid  in  this  we  augment  the  particle  equations 
with  the  corresponding  Vlasov  equation 


3  f  3  f  i  i<J>  -*  -i<K  3f 

W  +  P  3*  -  2  (ae  "  a  e  >3-7 
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(20) 


where  f(p,'l>,?)  is  the  distribution  function.  In  equilibrium,  f  =  fpCp). 
Linearizing  Eq.  (20),  assuming  exponential  dependences  of  the  first  order 
quantities  on  g,  and  using  Eq.  (15)  to  relate  the  field  amplitude  to  the 
perturbed  distribution  function  yields  the  dispersion  relation  for  v , 


!  =  Mil  1  I  dp  3f0(p) 


(v  +p  )  3  p 


(21) 


Here,  fg(p)  normalized  to  unity.  For  a  cold  beam  (fn  =■  ^(p-p^))  we 
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can  perform  the  integral  in  Eq .  (21)  to  obtain  1 

A 

For  large  beam  currents  I  »  1  we  anticipate  v  will  be  large  in  which 

case  we  can  use  the  large  argument  limit  of  E(v)  to  obtain 
2 

v(v+Pq)  =  -I/8tt  ,  which  is  identical  to  the  corresponding  one¬ 
dimensional  dispersion  relation.  Thus,  if  I  is  large  the  effects  of 
diffraction  are  negligible  mainly  because  fields  exponentiate  in  distance 


shorter  than  the  Rayleigh  length.  For  a  perfectly  matched  beam  (p 


0 


0) 


the  spatial  growth  rate  is  found  to  have  the  characteristic  third  root 

1/3 


dependence  on  current  v  =  ( — /  3/2)(I/8tt) 


If  I  is  small,  and  the  beam 


energy  is  perfectly  matched  we  expand  E(v)  for  small  to  v  to  obtain 
2  ,  _ _  -1 


(~I/87r)Jlnv  »  which  yields 

1/2 


-1/2 


'  8tt 


In  this  case  the  growth  rate  has  a  slightly  weaker  dependence  on  current 
A  1/2  A  1/3 

(between  I  and  I  )  owing  to  diffraction. 

A  detailed  analysis  of  the  dispersion  relation,  Eq.  (21)  including 
the  effects  of  beam  energy  mismatch  and  energy  spread  will  be  presented 
in  a  future  publication. 
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IV.  NONLINEAR  GUIDED  STATES 


We  now  consider  cases  where  the  field  amplitude  is  large  enough  to 
trap  particles  nonlinear ly  in  the  beat  wave  potential.  Further,  we 
assume  that  particles  execute  many  synchrotron  oscillations  before  the 
parameters  of  the  beat  wave  change.  To  treat  this  case  analytically  we 
let  the  field  be  expressed  in  terms  of  a  slowly  varying  amplitude  and 
rapidly  varying  phase  a(5  )  =  a^expCi  /^v(5') d?')  where  v (5 )  is  antici¬ 
pated  to  be  real  and  positive,  and  varies  slowly  with  5  in  response  to 
changes  in  the  wiggler  parameters.  In  the  particle  equations  we  define  a 

A  r 

new  phase  i|i  =  +  /  v(?')d£'  which  results  in  the  pair  of  equations 


9<J>  3  p  3 

=  p  v  ’ =  "3; 


—  aQcosi|>  , 


(22) 


where  we  have  taken  ag  to  be  real  and  we  have  neglected  the  deceleration 

a  to  lowest  order.  Equations  (22)  are  derivable  from  the  Hamiltonian  H 

2 

H  =  ^—  +  pv  +  agCosip  .  (23) 

Thus,  in  the  trapped  particle  regime  the  distribution  function  depends 
only  on  H,  f  =  f(H).  We  can  now  evaluate  the  source  term  in  Eq.  (14)  and 
obtain  the  expression  for  the  value  of  the  wave  amplitude  ag.  In  doing 
this  we  assume  that  v(5)  and  all  other  slowly  varying  quantities  are 
constant  on  a  length  scale  in  £  of  order  unity.  (In  fact,  for  small  v, 
because  of  the  slow  decay  of  kernal  we  must  make  the  more  stringent 
assumption  that  these  quantities  are  constant  on  a  scale  of  order  v  *.) 
We  then  perform  the  convolution  integral  in  Eq.  (14)  as  if  the  source  had 
a  simple  exp(iv^)  dependence.  Converting  the  p  integration  in  Eq.  (24) 
to  an  H  integration,  we  obtain 
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(24) 


*  [  2(H - ^  ~  aoC0^)] 

For  positive  values  of  ag  it  can  be  seen  from  the  Hamiltonian  that  deeply 
trapped  particles  will  have  cost|»  <  0.  Thus,  if  optically  guided  states 
exist,  the  average 


<cos^>  =  / 


#dH  f(H)cosij> 


2  A  1/2  » 

n[  2(H  -  —  -  agcos*)] 

must  be  negative.  Examining  E(v )  given  by  Eq.  (16a)  or  (16b)  It  is  seen 
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that  for  real  positive  v  E(v)  is  negative  and  thus  Eq .  (24)  is  consistent 
for  positive  Sq  and  steady  guided  states  exist.  The  existence  of  these 
guided  states  has  nothing  to  do  with  the  radial  profile  of  the  beam, 
rather  it  is  a  consequence  of  the  saturated  states  predicted  by  the  one¬ 
dimensional  model.  The  crucial  requirement  is  that  in  a  saturated  state 
the  phase  speed  of  the  current  perturbation  is  below  the  speed  of  light 
so  that  energy  can  not  be  radiated  laterally  from  the  beam. 

A 

Equation  (24)  indicates  that  even  for  very  weak  beams  I  «  1  guided 
states  exist.  This  is  due  to  the  logarithmic  divergence  of  E(v)  for 
small  v.  However,  the  filling  factor  for  these  states  will  be  very  small 
so  there  utility  may  be  limited. 

The  slow  evolution  of  the  guided  states  described  by  Eq .  (24)  can  be 
determined  by  returning  to  the  governing  equations,  and  restoring  the 
deceleration  a.  The  result  is  that  f  is  a  function  of  H  through  the 
conserved  particle  action  J  =  /  #p.  The  parameter  v  is  then  determined 
from  the  energy  conservation  relation 


.  1  2  3 E/3v  t  r  ... 

I<p>  +  r  a  — -  --if  d£  a  (£  ) 

U  EZ(v) 

where  the  average  of  <p>  is  found  to  be  simply  -v . 
we  can  write  the  energy  balance  equation, 


Thus,  using  Eq.  (24) 


1  3E  ,  “2  r 

2  <cos*|>  >  -  V  =  -  J  d(  o  . 

Recalling  the  definition  of  E  we  see  that  3E/3v  is  positive  and  decreases 

*  2 

with  v.  Thus,  assuming  the  dependence  of  <cosi|>>  on  v  to  be  weak,  we  see 
that  for  acceleration  (a  >  0)  v(5)  increases  with  £  corresponding  to 
improved  guiding  (i.e.,  a  larger  filling  factor).  For  deceleration, 

(a  <  0)  v(£)  decreases  with  £  corresponding  to  a  decrease  in  the  filling 
factor.  These  results  must  be  regarded  as  somewhat  qualitative  because 
of  the  assumption  of  adiabatic  change  in  the  wtggler  parameters.  A  more 
quantitative  evaluation  awaits  the  numerical  simulation  of  the  governing 
equations . 


V.  CONCLUSIONS 

In  this  paper  we  have  presented  a  model  system  of  equations  (the 
separable  beam  model)  which  is  capable  of  describing  a  wide  variety  of 
phenomena  associated  with  optical  guiding  in  free  electron  lasers.  The 
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model's  virtues  are  that  it  allows  for  an  exact  treatment  of  the  wave 
fields  and  collapses  the  description  of  the  FEL  interaction  back  to  a 
slightly  modified  version  of  the  much  studied  one-dimensional  model. 

Unfortunately,  the  key  step  in  obtaining  this  model  is  not  derivable 
from  first  principles.  That  is  the  separation  of  variables  in  the 
current  source  terra  is  not  rigorously  derived.  We  argue  however,  that 
the  details  of  the  radial  distribution  of  current  are  not  important  for 
beams  with  small  or  moderate  filling  factors.  These  are  beams  with 
insufficient  current  to  produce  exponentiation  of  fields  in  one  Rayleigh 
length.  In  some  respect  optical  guiding  in  these  beams  is  of  most 
interest  since  optical  guiding  in  strong  beams  is  more  assured. 

Finally,  we  suggest  that  due  to  the  exact  solvability  of  the 
radiation  field  in  the  separable  beam  model,  that  it  provides  a  test  by 
which  more  sophisticated  treatments  can  be  judged. 
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